HU-EP-05/42 



Laplacian modes probing gauge fields 

Falk Bruckmann" and Ernst-Michael Ilgenfritz 5 

Instituut-Lorentz for Theoretical Physics, University of Leiden, 
P.O.Box 9506, NL-2300 RA Leiden, The Netherlands 

b Institute of Physics, Humboldt University Berlin, 
Newtonstr. 15, D- 12489 Berlin, Germany 



We show that low-lying eigenmodes of the Laplace operator are suitable 
to represent properties of the underlying SU(2) lattice configurations. We 
study this for the case of finite temperature background fields, yet in the 
confinement phase. For calorons as classical solutions put on the lattice, 
the lowest mode localizes one of the constituent monopoles by a maximum 
and the other one by a minimum, respectively. We introduce adjustable 
phase boundary conditions in the time direction, under which the role of 
the monopoles in the mode localization is interchanged. Similar hopping 
phenomena are observed for thermalized configurations. We also investigate 
periodic and antiperiodic modes of the adjoint Laplacian for comparison. 

In the second part we introduce a new Fourier-like low-pass filter method. 
It provides link variables by truncating a sum involving the Laplacian eigen- 
modes. The filter not only reproduces classical structures, but also preserves 
the confining potential for thermalized ensembles. We give a first character- 
ization of the structures emerging from this procedure. 



1 Introduction 

The question what drives confinement and other nonperturbative phenomena of QCD at 
strong coupling is a long-standing one. In lattice gauge theory, the simulation of these 
effects at the observational level of correlation functions is well-established. However, 
to extract the relevant degrees of freedom of the QCD vacuum remains a controversial 
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problem which might not have a unique answer. The desire behind such attempts is to 
provide support for certain models like the instanton liquid or the dual Abelian Higgs 
model and to estimate their basic parameters. 

Abelian and center projection have been used to focus on objects like Abelian mag- 
netic monopoles and center vortices, respectively, which are thought of as localizing 
certain embedded solutions that, strictly speaking, would exist only in the presence of 
Higgs fields. Despite the ambiguity in their definitions, these degrees of freedom are 
used especially for confinement-related scenarios of the QCD vacuum. Instantons as 
selfdual solutions, on the other hand, can be conveniently related to chiral symmetry 
breaking. Field excitations resembling instantons have indeed been observed as the re- 
sult of cooling. It has been objected that these methods modify the field configurations 
in an uncontrolled way such that the observed excitations could actually be fake and 
do not represent the relevant nonperturbative fields. Other smoothing techniques also 
reveal lumps of action which are usually interpreted as instantons ^ Hj • In order to 
relate these structures to confinement, some additional degrees of freedom seemed to be 
necessary. 

Fermionic modes capture chiral and topological aspects of lattice gauge theory pro- 
vided a Dirac operator with good chiral properties is implemented. Via the index the- 
orem the number of (left-handed minus right-handed) zero modes gives the total topo- 
logical charge. Using a definition of the topological charge density based on the overlap 
Dirac operator [3], non-zero modes are contributing to the topological charge distribu- 
tion, too. Analyzing the latter, global coherent structures of lower dimension have been 
identified @]. 

It has become customary to use the localization of the fermionic modes as a means 
to probe the vacuum structure. A specific scaling law of the inverse participation ratio 
(IPR, see below) of the low-lying fermionic modes w.r.t. lattice spacing and volume seems 
to be able to recognize the codimension of the underlying gauge field structure 
However, to interpret the result, the latter has to be subject to a model, for instance 
one of those mentioned above. 

The fact, that fermion zero modes are localized to classical objects in smooth back- 
grounds, has lead to a thorough investigation of the local properties of zero and near-zero 
modes also in equilibrium backgrounds. Their small 'energy' eigenvalue should suppress 
contributions from large momenta and hence these modes are less sensitive to UV fluc- 
tuations. Very recently such a programme has been carried out using zero modes in the 
adjoint representation 7 . 

The trick of modifying the boundary conditions for the fundamental fermions by a 
complex phase has been introduced as a general tool in |Hj guided by the knowledge of 
calorons. These are instantons at finite temperature, i.e. on i? 3 x S 1 . Calorons have 
become attractive over the recent years because - when taken with nontrivial holonomy 
(see below) - they account for a Polyakov loop not in the center of the gauge group, 
as is the case on average in the confined phase. Furthermore, calorons contain (gauge 
independent) magnetic monopoles, see Fig. ^(a), which realize the scenario of fractional 
charge objects (also called instanton quarks JJ). For recent progress on calorons, both 
in the continuum and on the lattice, see |12j . 
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Figure 1: (a) The action and topological density (not distinguishable, both multiplied 
by 400) and the Polyakov loop shown along the line connecting the constituents 
in a large caloron on a 16 3 • 4 lattice, (b) The modulus of periodic (boldface 
line) and antiperiodic (thin line) fermion zero modes for the same caloron. 

The caloron zero mode with different boundary conditions is correlated to different 
constituent monopoles ^3], see Fig. ^(b) 1 . In a similar way the zero modes on equilib- 
rium configurations have been observed to localize to different locations on the lattice 
when scanning through the boundary conditions [B]- This effect has been reported even 
for symmetric lattices representing 'zero temperature' |14j . Given the intuition from 
calorons, one expects these modes to detect carriers of topological charge, including 
such of fractional charge. However, to support this one would need additional evidence 
that the underlying structures actually have fractional or integer charge. A mechanism 
analogous to Anderson localization in a random potential |15j has been proposed as an 
alternative explanation for the localization and hopping of the fermionic modes. 

In this paper we are going to study the analytic power of the low-lying eigenmodes of 
the gauge covariant Laplace operator. They have been first discussed with the aim to fix 
a particular gauge, the Laplacian gauge [TSJ. Recently they have been investigated with 
respect to their localization behaviour [T2 a to shed additional light on the QCD vacuum. 

Our interest concentrates on the local behaviour of these modes testing three ideas. 
The first idea is whether the Laplacian modes are sensitive to the constituents of the 
caloron. Unlike fermion zero modes, the Laplacian modes are not incorporated in the 
ADHM-Nahm formalism 2 that describes the caloron solutions analytically. So we inves- 
tigate them numerically on the lattice, on which caloron configurations can be put with 
very good control. We find indeed that the Laplacian modes can detect the monopole 
constituents by their extrema. Furthermore, we enforce adjustable phase boundary con- 
ditions on the Laplacian modes as a function of which they are hopping in a similar way 
as the fermion modes do. 

Most of our studies are concerned with eigenmodes of the Laplacian in the fundamental 



1 We thank Dirk Peschka for providing the fermion zero mode for the (numerical) caloron. 

2 although the Greens function is 
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representation. For comparison, we also explore adjoint modes, including those with 
antiperiodic boundary conditions. 

At this point we want to emphasize that the Laplacian modes are advantageous com- 
pared to the fermionic modes in that they have neither chirality nor doubler problems. 
Hence, the straightforward translation of the continuum Laplace operator to the lattice 
can be used for the purposes we have in mind. We will mostly study the lowest 'en- 
ergy' eigenmode which has no topological origin (and therefore is not expected to give 
information about the topological charge). 

The other idea is how Laplacian modes could be used to reflect properties of equilib- 
rium configurations. As we will demonstrate, they do this in the 'conventional way' by 
being pinned to some points on the lattice and hopping between several such locations 
under a change of the boundary conditions. 

The relation between the localization mechanisms for classical and thermalized back- 
grounds is not straightforward. As we will show, the lowest Laplacian eigenmode for the 
caloron is similar to a modified wave (in contrast to the exponentially localized fermion 
zero mode). Its minimum, but even more the maximum for small calorons, is not a 
perfect tool, whereas in the thermalized backgrounds the maxima obviously determine 
the localization and divide the set of boundary conditions into intervals. 

The third main aspect of the present work is the introduction of a general and hope- 
fully powerful method to apply a low-pass filter based on the Laplacian eigenmodes to 
generic equilibrium configurations. We were inspired by Gattringer's earlier work who 
has constructed a smoothed field strength tensor based on fermionic modes |18j . Our 
procedure is not limited to a particular observable, but aims to reconstruct the link vari- 
ables filtering out UV fluctuations. It relies on the idea of truncating a sum involving 
the Laplacian ('harmonic') modes which would give back the original links exactly. The 
only free parameter (besides the number of modes used in the truncation) is the phase 
angle in the boundary condition. We will show that in order to optimize the low-pass 
filter for our circumstances (confined phase, nontrivial holonomy calorons), this angle 
needs to be chosen halfway between periodic and antiperiodic boundary conditions. 

For the filtered configuration we measure the Polyakov loop, the topological charge 
and action density. Again, we present the calorons as a testing ground and find, indeed, 
that the selfdual monopole constituents are reproduced. Qualitative agreement with the 
original configuration is found already with a surprisingly small number of modes. 

What is even more interesting is that this method approximately preserves the string 
tension when applied to Monte Carlo configurations representing some finite temperature 
below the deconfinement phase transition. Apparently, the Laplacian modes capture 
enough of the long range disorder. This fact justifies to study the filter in more detail, 
especially the emerging tomography of the configurations in equilibrium. The structure 
of the action density of the filtered configurations includes narrow peaks, for which we 
do not have a final interpretation yet. They might be close to the gauge singularities 
found in the Fourier-filtered Landau gauge fields |19j . Despite the appearance of the 
peaks we will point out a similarity to smearing/cooling in an early stage. 

The paper is organised as follows. In the next section we will give the definition 
as well as some properties of the lattice covariant Laplacian and briefly summarize 
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the knowledge about calorons. Then we investigate low-lying Laplacian eigenmodes in 
caloron backgrounds. In Sect. 0] we show the typical hopping of the Laplacian eigenmodes 
in the background of thermalized configurations at finite temperature. The new filter 
method and its features are discussed in Sect. |SJ We end with some discussion and an 
outlook. We will stick to the gauge group SU (2) throughout this paper. 



2 Preliminaries 



2.1 Definition and properties of the Laplace operator 

We consider the gauge covariant Laplace operator 

D 

A t = E [ U ?( x ) s *+H,v + u t ab (y) S *-fi,y - 2Sa * S *v\ , a,b=l,2, (1) 

in the background of a given configuration of lattice links U^(x) in D = 4 dimensions and 
in the fundamental representation. It is a hermitean (and non-positive) matrix of size 
N x M where N is the number of lattice sites times the dimension of the representation 
(here two). 

We use the ARPACK package fJU] to solve for up to 200 (out of M > 10 5 ) eigenvalues 
and eigenmodes in 

- A ^</>n,c(y) = Vc0n,c( X ) 110 SUm 111 U aIld ( 2 ) 

We allow for complex phase boundary conditions in the time-like direction 

(t> c {x A + N A ) = e 2 -^ c (x 4 ). (3) 
A way to implement this is by writing 

<pdx)=e- 2 * iX4</N4 <t>d*), ( 4 ) 

which is now fully periodic but solves the Laplace equation © with the replacement 
U4 — > Uie 2 ^^ 4 - in A. Effectively, this promotes the link to an element of U(2), but 
does not change any contractible Wilson loop. Correspondingly, the gauge field receives a 
constant identity component, which does not modify the field strength. In the numerical 
computations we use (p and transform back to <j) by virtue of 

The spectrum of the Laplace operator is subject to two symmetries. The charge 
conjugation (which is special for SU{2)) 

<_ c = e°%, A^_ c = A n>c (5) 

relates two modes with the same eigenvalue but opposite boundary conditions, as indi- 
cated by the index — £. In particular, the spectrum is two-fold degenerate for periodic 
(£ = 0) and antiperiodic (£ = 1/2) boundary conditions. Notice that 4>' is automati- 
cally orthogonal to <fi. We will therefore restrict ourselves to £ £ [0, 1/2] without loss of 
generality. 
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The authors of Ref. have found a relation between the lower and upper end of the 
spectrum, 

^-n,( = (-l) E ^"^(. ^_ niC = 4L>-A n , c . (6) 

We will refer to this symmetry as the staggered symmetry. It is only for even numbers 
Nn of lattice points in all directions that it (obviously) preserves both the boundary 
condition Eq. @ and the periodicity in the space-like directions. Thus, even numbers 
N/j, will be used throughout. This symmetry, since it flips the sign of 4> at every other 
point, is restricted to the discrete lattice (otherwise the spectrum of the continuum 
Laplacian would be bounded from above as well as from below). 

2.2 Laplacian modes in vacuum backgrounds 

In order to illustrate the interplay of the boundary condition angle C an d the Polyakov 
loop 

N 4 

V{x) = J] U 4 {x,x 4 ) (7) 

X4=l 

with the spectrum of Laplacian modes, we now discuss the latter for vacuum configura- 
tions. We choose all time-like links to be identical and to belong to the Abelian subgroup 
consisting of diagonal matrices U±(x) = exp(27ria<T3/iV4) which leads to a constant but 
adjustable Polyakov loop 

-tvV(x) = cos(2vtq) , (8) 

whereas the space-like links are set equal to the identity. 

The Laplace equation in these backgrounds can be solved by considering the upper and 
lower component separately and a simple product ansatz <f>{x) = JT exp(27ri/3 /t x M ). The 
periodicity requirements give = nj/iVj, = {n^ + Q) /N A with integers £ [0, N^ — l]. 
The eigenvalues are 



,n 4 C±a y-^ rii 

, + > 1 — COS27T 



(9) 



One can immediately read off that the dependence on £ is trigonometric and shifted by 
the Polyakov loop parameter a, while the spatial part shifts the spectrum and contributes 
to the degeneracy. 

In Fig.|2]we give the vacuum spectra for a = and a = 1/4, which will be of particular 
interest below. We plot the 38 lowest modes 3 and indicate the degeneracy of the bands. 
They are reflected at the boundary of the 'Brillouin zone' (where the eigenmodes are at 
least two-fold degenerate as they should) and then cross each other. 

Aspects of these spectra, especially the symmetry around £ = 1/4 for a = 1/4, will 
be reproduced by Laplacian modes in calorons and, to some extent, in thermalized 
backgrounds. 

3 in order to completely fill the three lowest-lying bands at £ = 
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Figure 2: Spectrum of the Laplacian in vacuum backgrounds with (a) traceless (a = 1/4) 
and (b) trivial Polyakov loop (a = 0) on a 16 3 • 4 lattice. Plotted are the 38 
lowest eigenvalues as a function of the angle £ in the boundary condition. The 
numbers indicate the degeneracy of the bands. The unit on the vertical scale 
is 0.152 = 2(1 - cos(2vr/16)), cf. Eq. ©. 



2.3 Calorons - in the continuum and on the lattice 

Calorons are (anti-) self dual instanton solutions on R 3 x S . We will be particularly 
interested in caloron solutions with the so-called holonomy being maximally nontrivial jUJ 
[TIT] . That is tr Poo/2 = cos(27raoo) = 0, where Voo is the Polyakov loop at spatial infinity. 
These, rather than the old Harrington-Shephard solutions |U (with Voo = ±12) > shall 
be of relevance for the confined phase, where the trace of the Polyakov loop vanishes on 
average. The authors of j^2] have computed the contribution of calorons to the effective 
potential driving the Polyakov loop to that value. 

The nontrivial holonomy gives rise to a symmetry breaking SU(2) — * U(l). There- 
fore, it is plausible that calorons can be described by a monopole and an antimonopole 
(calorons of charge k consist of |fc| monopoles and \k\ antimonopoles). Their masses are 
given by the eigenvalues of the holonomy and equal for maximally nontrivial holonomy. 

The moduli space of these solutions contains both small and large calorons, where 
the size is compared to the extension of the compact direction. Large calorons have 
two almost static action density lumps, which are identical for our choice of holonomy. 
They merge for small calorons and this results in a strong time-dependence of the action 
density (just like for conventional instantons). 

The (untraced) Polyakov loop acts like a Higgs field (exponentiated to give an element 
of the gauge group). It passes through I2 and — 12 near the core of the monopoles, where 
the symmetry is restored, see Fig. ^ (a). This dipole persists even for small calorons 
within the single action density lump. 

The monopoles make up the topological charge with the help of the so-called Taubes 
winding [23]: one of the monopoles performs a full rotation in the unbroken U(l)- 
subgroup relative to the other monopole when completing a full period in the time-like 
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direction. For a gauge invariant statement one has to connect the field strength at the 
different monopole cores by a Schwinger line [21] • In the periodic gauge that we use 
for the calorons put on the lattice, the link variables U^x) are static at the V = I2 
monopole, while they rotate around the holonomy direction T3 at the V = — 12 monopole. 

The fermion zero mode in the caloron background follows the action density in that 
it has a maximum at one constituent monopole ^3] plus a zero near the other one 
|25| . The first fact can be understood from the Camas index theorem [2111, which also 
explains why the zero mode hops with the angle £, namely to the other monopole when 
the boundary condition is changed from periodic to antiperiodic. When the phase £ 
equals the holonomy parameter aoo, the zero mode 'sees' both monopoles, decaying 
with a power-law instead of exponentially. The zero of the zero mode is connected to 
the nontrivial caloron topology. 

On the lattice, calorons were first obtained by cooling in Ref. |24j . where the typi- 
cal behaviour of the action density, the Polyakov loop and the fermion zero mode was 
confirmed. Furthermore, the Polyakov loop averaged over the low-action region of the 
lattice was proposed to play the role of the asymptotic holonomy in the infinite contin- 
uum. As for generic equilibrium configurations, recent smearing studies revealed clusters 
of topological charge, that - according to their content of Abelian monopoles and the 
corresponding behaviour of the Polyakov loop - have charges close to either ±1 or ±1/2, 
resembling calorons and their constituents, respectively [2*7] . 

An alternative possibility to analyze calorons on the lattice is to evaluate the discrete 
parallel transporters from the continuum gauge field and to cool the emerging lattice 
configuration by a few steps in order to adapt it to periodic boundary conditions in space. 
This was first done in [2Hj for the case of gauge group SU (3). Because of lattice artefacts, 
the use of improved or even overimproved cooling |29| is advantageous when aiming at 
large calorons. We will mainly use this approach to generate caloron backgrounds since 
it permits to control the locations of the constituent monopoles. 

3 Laplacian modes in caloron backgrounds 
3.1 Eigenvalue spectra 

The lowest-lying Laplacian eigenvalues in a background of a small and a large caloron, 
both with maximally nontrivial holonomy, are displayed in Fig. 01 (a) and (b), respec- 
tively. By construction, the trace of the Polyakov loop of these lattice configurations 
vanishes outside of the action density lumps of the constituent monopoles. Accordingly, 
the spectra of the Laplacian are similar to those of a vacuum background with the cor- 
responding Polyakov loop with a = 1/4 in Eq. (jHJ (see Fig. [21(a)). However, the zero in 
the vacuum spectrum at £ = 1/4 is lifted, the degenerate bands are splitted and some of 
the level crossings are inevitably avoided. Such effects are known form quantum systems 
at weak coupling (with the Laplacian playing the role of a Schrodinger operator). In 
this respect, the caloron backgrounds act like a small perturbation. 

In order to characterise the localization of the eigenmodes by a global quantity, we 
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Figure 3: Dependence of the 38 lowest-lying eigenvalues on the boundary condition for 
(a) a small caloron and (b) a large caloron, both with maximally nontrivial 
holonomy on a 16 3 • 4 lattice. 



will use the inverse participation ratio (IPR) 

I(<f>) = Vol- J> 2 (x), p(x) = |0| 2 (x), Vol = NxN 2 N 3 N 4 . (10) 

X 

A constant profile has the minimal 1 = 1, whereas large IPR's signal strong localization 
(with the lattice (5-function saturating the upper bound / = Vol). For the IPR's of 
fermionic modes in caloron backgrounds we refer the reader to Ref. |3()j . 

In Fig. |1] we zoom in to a region of (almost) crossing eigenvalues. In contrast to the 
vacuum spectrum, the eigenvalues of the second and third state are repelled from each 




Figure 4: Zooming in to (a) the eigenvalues number 2 to 8 and (b) the corresponding 
inverse participation ratios for the small caloron shown in the previous Fig. 
01 (a). For comparison also the IPR of the lowest mode (closest to / = 1) is 
shown. 
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Figure 5: Low-lying spectra (lowest 15 modes) for charge 2 calorons, (a) the case of 
pairwise well-separated constituents, (b) the case of like-charge constituent 
monopoles merged forming double-monopoles (rings, see the text and Ref. 

m). 



other. As a remnant of the crossing these modes exchange their IPR's around that £, 
which indicates that in the next (^-region each mode is similar to the complementary mode 
in the previous region. Outside of the crossing regions the IPR's are hardly changing. 

At some other crossing points the eigenvalues come very close and the switch to the 
other eigenvalue branch is an instantaneous one even within the enlarged resolution in 

The drastic changes in the IPR's illustrate this again. One might speculate whether 
some of the level crossings are exact ones in the continuum limit. They may be governed 
by the quantum number corresponding to the axial symmetry of the caloron. However, 
as the spatial volume increases, Ni — > oo, many bands will approach each other to form a 
continuous spectrum and to investigate the existence of discrete eigenstate would require 
a more detailed study. 

In Fig. [3 the spectra of two extreme cases of charge-2 calorons (32] are shown. The 
four constituent monopoles are maximally separated in the first example, whereas in the 
second case the like-charge monopoles sit on top of each other forming rings (see Fig. 7 
and Fig. 5 of Ref. [SI]). 

In an overall view, the spectra are similar to those of charge-1 calorons. However, 
the eigenvalues are shifted upwards, especially for the second example (Fig. |S] (b)). The 
band structure is different as well, with three close eigenvalues at the bottom of that 
spectrum. The other case (of maximally separated constituents in the charge-2 caloron) 
seems to have curves with different curvature (seen in the middle of Fig. |S] (a)) which 
upon a closer look turn out to have their minima slightly away from £ = 1/4. Although 
these details are interesting, it is not clear how one can infer, for instance, the topological 
charge of the background configuration from these features. 

In all the gauge field backgrounds studied so far, the IPR of the lowest eigenmode is 
close to / = 1. It means that this mode is quite spread out. For the excited states the 
IPR takes values up to / = 3.5. 
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3.2 Eigenmode profiles 

In this subsection we study the local properties of the Laplacian modes for a caloron 
background and show how they can reveal the underlying monopole constituents. 

The modulus of the lowest mode for a large caloron is plotted in Fig.^fbr different 
boundary conditions. This figure should be compared to Fig. (a) (for the background 
field) and Fig. U(b) (for the fermion zero modes). One can see that the lowest Laplacian 
eigenmode with periodic boundary conditions has a maximum (located at X3 = 14) 
near the monopole with a positive Polyakov loop (at x% = 12). Furthermore, it has a 
minimum (at X3 = 6) near the monopole with negative Polyakov loop (at X3 = 4). Thus, 
it approximately localizes the two constituents by virtue of a minimum and a maximum. 
However, both are not very pronounced and occur with a shift of up to 2 lattice spacings 
(compared to a time-extent of N4 = 4 lattice spacings) relative to the locations of the 
action density lumps. 

We have also investigated excited eigenstates of the Laplacian and found that the first 
excited one has a higher maximum, but both minimum and maximum are now shifted 
in the opposite direction. Even higher states seem to be rather sensitive to the finite 
spatial volume. 

What can also be read off from Fig. is that the minimum and the maximum of the 
lowest mode move to the other constituent upon changing the boundary condition from pe- 
riodic to antiperiodic. This symmetry can be understood in the continuum. The caloron 
action density is invariant under an antiperiodic gauge transformation that exchanges 
the locations and masses of the constituent monopoles and inverts the holonomy. The 
latter two replacements are uneffective for a caloron at maximally nontrivial holonomy. 
Hence, for the Laplacian mode this gauge transformation reflects the profile function at 
the caloron's center of mass (around £3 = 8). Furthermore, it replaces £ by £ + 1/2 (due 
to the antiperiodicity) which by charge conjugation is equivalent to 1/2 — Q. Thus the 
periodic mode turns into the antiperiodic one. This particular caloron symmetry also 
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Figure 6: The modulus |0| of the lowest Laplacian mode with periodic (boldface line), 
antiperiodic (thin line) and intermediate (dashed line) boundary conditions 
along the line connecting the constituents of a large caloron on 16 3 • 4, to be 
compared with Fig. ^ 
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Figure 7: Space-space plots of the modulus of the lowest eigenstate for a large caloron 
on a 32 3 • 4 lattic, with (a) periodic and (b) intermediate boundary condition. 



explains the mirror symmetry of the spectrum in Fig. [3] at £ = 1 /4. 

In this respect the lowest Laplacian eigenmode resembles the fermion zero mode. 
However, the interpolation between the two extreme boundary conditions is different: 
The minimum moves through the center of mass, whereas the maximum decreases and 
goes through the 'boundary' at 'infinity'. To demonstrate this we have included the 
intermediate boundary condition £ = 0.25 in the figure. 

In order to clarify how finite volume effects influence these findings, we have doubled 
both the spatial extension of the lattice and the size of the caloron, which now has its 
constituents at x\ = X2 = 16, X3 = 8 and 25. The lowest eigenstate with periodic 
boundary condition, shown in Fig. [3(a), has a slightly more pronounced maximum, still 
shifted, to X3 = 26. Around the maximum the modulus is actually close to spherically 
symmetric as expected since the constituents are quite well separated compared to the 
typical scale N4 = 4. The minimum is shifted as well, namely to £3 = 10. 

This configuration also clarifies the behaviour of the minimum under the change of 
boundary conditions. Fig. (b) shows that for an intermediate boundary condition 
the minimum extends to a valley between the minimum locations corresponding to the 
periodic and antiperiodic modes, respectively. 

It is instructive to see how the lowest Laplacian eigenmode looks like in the background 
of a smaller caloron. The configuration that we employ for this demonstration has only 
one lump of action density, but the typical dipole behaviour of the Polyakov loop, the 
latter being close to — 12 at X3 = 15 and close to I2 at X3 = 18. The modulus of the 
lowest Laplacian eigenmode has a strong gradient at the location of the lump, too, see 
Fig- El (a)- As a matter of fact, the minimum (2:3 = 15) reflects the Polyakov loop 
minimum quite well, while the maximum (X3 = 23) is shifted far outwards. 

The maximum is also less pronounced compared to the one of the large caloron. 
Actually all the presented profiles of the lowest Laplacian eigenmodes away from the 
monopoles become close to the value an entirely constant mode would have on the 
corresponding lattice (0.0078 for 16 3 • 4, 0.0028 for 32 3 • 4). Therefore, the Laplacian 
modes should be viewed as a wave with local modifications rather than a localized state. 
The rationale behind that might be the absence of an effective mass that would localize 
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Figure 8: (a) Space-space and (b) space-time plot of the modulus of the lowest eigenstate 
with periodic boundary conditions on the background of a small caloron on a 
32 3 • 4 lattice. 



the modes exponentially. For the fermion zero modes this role is played by (the difference 
of the boundary condition and) the eigenvalues of the holonomy. 

Moreover, the maximum of the Laplacian mode is almost static, even for the small 
'instanton-like' caloron, as can be seen in Fig. [H] (b). The minimum, on the other 
hand, shows up in a particular time-slice, exactly where the action density is maximal. 
This is another sign that, from a practical point of view, the use of the minimum is 
superior compared to the maximum when looking for small calorons (or instantons). 
The time dependence of the minimum becomes much weaker for large calorons as we 
have observed (but not shown here), reflecting the more static character of the solution. 
In the continuum, Laplacian eigenmodes in nontrivial backgrounds will have a zero of 
topological origin, like fermions do |25j . 

To summarize, the assignment of a maximum and a minimum in the modulus of the 
lowest Laplacian eigenmode to the constituents in a caloron is similar to the behaviour 
of the fermion zero mode. Yet the profile is less localized and shifted. For small calorons 
it is recommendable to use the minimum, as it was done in Ref. [33], to find instantons 
on the lattice (in the continuum, the coincidence of the minimum with the instanton 
center was shown analytically in [31] )• The shift observed in that reference nicely agrees 
with our findings. 

By inspecting the components of the Laplacian modes one finds signatures of the 
Taubes winding. We plot the real and imaginary part of <P^Zq' 2 in Fig. [2] showing the 
time dependence at both monopoles. At one monopole core the components are fully 
static, Fig. 03(b). At the other monopole, carrying the Taubes winding, they oscillate 
over one period, however, they seem to contain an additional constant part, Fig. [5] (a). 

The general picture of Laplacian modes described so far passes over to the caloron 
examples of charge Q = 2. The periodic and antiperiodic modes possess maxima and 
minima near the corresponding constituents. However, the ring structure in the action 
density of two overlapping like-charge monopoles is not resolved by the Laplacian modes, 
in contrast to the fermion zero modes. 
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Figure 9: Components of the lowest eigenstate with periodic boundary conditions, in 
both the fundamental (full lines, real and imaginary part of both components) 
and the adjoint representation (dashed line, the T3-component as dotted line), 
as function of £4 at (a) the point x\ = xi = 16, X3 = 8 showing the Taubes 
winding and (b) the point x\ = X2 = 16, X3 = 25 (no winding). Two pairs 
of lines fall on top of each other in plot (b): real and imaginary part of the 
second fundamental component and the first two adjoint components. 



3.3 Adjoint representation 

In this subsection we discuss (in short) the properties of the lowest Laplacian eigenmode 
in the adjoint representation. In the Laplacian, the fundamental links U® (x) are simply 
replaced by the adjoint ones ti((j A U ll (x)a B Uj l (x)) /2 with indices A, B running from 
1 to 3. The eigenfunctions <ft A (x) are rea l i n this representation and so there is no 
continuous phase available to modify the boundary conditions as this was possible in 
the fundamental representation, Eq. (J3J) and (JIJ). Nevertheless, we have included the 
possibility of antiperiodic boundary conditions in the computations. 

In caloron backgrounds we basically reproduce the findings of [S3]) namely that the 
constituent monopoles are localized by minima in the modulus of the lowest-lying ad- 
joint Laplacian eigenmode with periodic boundary condition. In other words, the single 
adjoint mode 'sees' monopoles of both kinds simultaneously. As Fig. EH shows, these 
minima have no shift problems and start to join for a small caloron. That the lowest ad- 
joint Laplacian eigenmode on an instanton background has a zero of second order at the 
instanton location has been shown in the continuum in [34| . Zeroes in this mode underly 
the construction of Abelian projected monopoles in the Laplacian Abelian gauge |36| . 
Thus, for calorons these monopoles coincide with the (gauge independent) constituent 
monopoles. 

The lowest adjoint Laplacian eigenmode with antiperiodic boundary condition devel- 
ops an extended static zero sheet between the monopoles (similar to the lowest funda- 
mental eigenmode with intermediate boundary conditions). The corresponding profile 
is shown in Fig. a s a dashed line. 

This feature can be of interest for detection purposes, too. In order to illustrate this, we 
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Figure 10: The lowest adjoint eigenmodes (with periodic boundary condition) for the 
large and the small caloron on a 32 3 • 4 lattice. The dashed line shows the 
lowest adjoint eigenmode with antiperiodic boundary condition for the large 
caloron. 

have investigated a semiclassical case. The gauge field configuration has been obtained 
by cooling down to a plateau with an action of 1.92 instanton units and a topological 
charge of —1 As has been described in that reference, there are one selfdual and 
three antiselfdual monopoles, distinguished by color in Fig. ^J(a). The Polyakov loop 
at one of the cores is found close to I2 and close to — 12 at the three others. Isosurfaces 




xl xl 



(a) (b) 

Figure 11: The semiclassical configuration described in the text: (a) isosurfaces of the 
topological density (positive as dark green, negative as light red), (b) isosur- 
faces locating the minima of the modulus of the lowest adjoint modes, both 
for periodic (pointlike, balls corresponding to the constituents seen in (a)) 
and antiperiodic boundary conditions (one-dimensional network). 
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(at small value) of the modulus of the lowest adjoint Laplacian eigenmode with either 
boundary condition are shown in Fig. ^J(b). As for the simple (dissociated) caloron, 
the periodic mode has static minima at the monopole locations, whereas the regions of 
small modulus of the lowest antiperiodic mode form a network connecting them. 

The lowest adjoint Laplacian eigenmode with periodic boundary condition actually 
reveals the Taubes winding in a very clear manner, cf. Fig. |5J At the 'rotating' monopole 
(a) the t\ and T2 components rotate around the holonomy subgroup generated by T3, 
whereas at the other monopole (b) all components are static. For the antiperiodic 
lowest mode all components perform half a rotation everywhere (to account for the 
antiperiodicity) with the T3-component being suppressed. 

4 Hopping of the lowest Laplacian mode for thermalized 
configurations 

In the last section we have shown that the lowest Laplacian eigenmodes (fundamental and 
adjoint) reflect certain properties of smooth classical gauge field backgrounds. Now we 
explore these modes as an analyzing tool for thermalized gauge field configurations and 
concentrate on how the dependence on boundary conditions gives additional information. 

As a set of thermalized background configurations we take an ensemble of 50 configu- 
rations on a 16 3 -4 lattice, generated by Monte Carlo heat bath sampling with the Wilson 
action at (3 = 2.2. It represents the confining phase at finite temperature, namely at 
T ~ 0.75 T c (deduced from a(0)a 2 = 0.22 for our f3 jSH and T c /y^(&j = 0.709 for SU{2) 
YM theory 0g|). 

Fig. El (a) shows a typical spectral flow of the 15 lowest-lying Laplacian modes with 
the boundary condition angle £. The first outstanding feature to notice is that the 
eigenvalues themselves are much bigger 4 than for the smooth backgrounds considered so 
far. This is most naturally ascribed to the ultraviolet noise present in the background. 
The latter has also removed any remnants of the vacuum (or caloron) band structure in 
this plot. Still, the typical near-crossing points are present, where again the IPR signals 
big rearrangements in Fig. 1121 (b). see e.g. the behaviour of the second and third mode 
around Q = 0.12. 

Inspecting several independent configurations in the confined phase, the lowest eigen- 
value takes on its smallest value around £ = 1/4 (in this respect the configuration in Fig. 
12 is not a typical one). This reflects the fact that the average Polyakov loop is close 
to traceless. On the contrary, in the deconfined phase all low-lying eigenvalues become 
minimal very close to £ = (and are much bigger at £ = 1/2) with the spectrum grossly 
similar to Fig. [21(b), in accordance with the asymmetric Polyakov loop distribution in 
that phase. 

The main observation in the context of Laplacian modes for thermalized backgrounds 
is the effect of 'hopping' of these modes with changing £. The IPR of the lowest mode 
already gives a hint on the existence of what we will call the '^-intervals'. In the example 



4 With our lattice spacing the eigenvalue unit in Fig. 1121 fa) is 1/a 2 = 0.87 GeV 2 . 
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Figure 12: Laplacian modes for a configuration in the confined phase. Plotted over the 
boundary condition £ are (a) the 15 lowest-lying eigenvalues in lattice units, 
(b) the IPR of the 6 lowest-lying modes (with corresponding symbols); (c) 
the value of the modulus of the lowest eigenmode at its global maximum and 
minimum (multiplied by a factor 100), (d) coordinates of the global maximum 
of that mode, displaying the typical (-intervals. 



of Fig. El (b) there are three such intervals, inbetween which the IPR has a minimum. 
It means that the lowest mode delocalizes there in order to rearrange itself. This is 
confirmed by the inspection of the global maximum in Fig. EES (c): the value of the 
modulus at the maximum is minimal at the two transition points ( = 0.300 and ( = 
0.415. The main feature of the ("-intervals is the pinning down of the lowest mode to 
particular locations. That means that the coordinates of the global maximum have 
constant values within the (-intervals and jump at the transition points, as it is clearly 
visible in Fig. El (d). 

For some configurations we have observed two distinct minima in the lowest eigen- 
value as a function of (, which turned out to be another signal for the existence of the 
(-intervals. 
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In the configuration discussed above the global maximum jumps over spatial distances 
of 8.5 and 12.7. We will restrict ourselves for the moment to a spatial analysis, since 
the Laplacian modes are almost static (see below and Fig. 116(1 . We have inspected 
50 configurations with 31 boundary conditions and recorded the spatial jumps. Most of 
them (160) are smaller than two lattice spacings, which can be viewed as a discretization 
error, where coordinates of the maximum change by one lattice spacing. We have not 
considered these minimal jumps further. 

There remained 120 jumps over at least two lattice spacings. Their distribution is 
plotted in Fig. 1131 It is still dominated by small jumps, however, the jump distance also 
reaches values almost as high as the maximally possible distance 13.9 = (16/2 • 
There seems to be no correlation of this distance to the £-value at which the jump 
occurs. Instead, the fact that the jump distribution is rather flat around half the linear 
extension of the lattice (here 8), may point to a random distribution of pinning centres, 
as suggested for the hopping of fermions in ^3] . In order to make such a statement more 
precise, one would need a model also for the small jumps. 

Having underlying calorons in mind, small jumps would be related to small calorons 
(and larger calorons might actually be suppressed in the same way as large instantons). 
However, calorons of holonomy close to maximally nontrivial (as expected for the con- 
fined phase) would prefer one jump at £ ~ 1/4. On the other hand, out of the 50 
thermalized configurations we find only 11 configurations with only one jump, but 16 
with two, 15 with three and 8 with even four jumps, respectively, and the jumps are 
scattered between ( = and ( = 1/2. 

In the last section we have argued that actually the minima of the modulus of the 
lowest Laplacian mode are good markers for classical topological objects. In the thermal- 
ized background, however, the global minimum is not as stable as the global maximum. 
Fig. ^1 (c) shows that for our example its value is changing smoothly only in a certain 
fraction of the first ^-interval. For most values of £ the value is highly fluctuating and 
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Figure 13: Statistics of the spatial distance of jumps from 50 thermalized configurations 
(160 jumps occured with spatial distance between 1 and 2 lattice spacings 
but were not considered further). 
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so is the location (not shown) . Furthermore, the number of local minima is typically an 
order of magnitude bigger than the number of local maxima. Therefore, it is difficult to 
employ those minima to eventually localize background instantons or calorons within a 
thermalized gauge field. 

In Fig. El we inspect the lowest mode locally. As can be seen immediately, the 
modulus of the lowest mode is quite smooth which justifies the expectation that these 
modes do not possess UV fluctuations. Horizontally in the figure we have plotted three 
lattice planes that contain the locations of the global maximum taken in the different 
£ intervals. The precise (- values representing the intervals are chosen by the demand 
for a maximal IPR. Following the maxima one can see that they either remain as local 
maxima or disappear in the other (-intervals. In other words, the hopping effect comes 
about by local maxima turning into global ones in particular ("-intervals, as visualized 
in Fig. El 

For comparison we show in the first and fifth row of Fig. El the lowest adjoint mode 
with periodic and antiperiodic boundary condition, respectively. These modes are corre- 
lated to the fundamental ones with the closest boundary conditions. Hence the adjoint 
Laplacian modes are also subject to hopping (when one allows for antiperiodic boundary 
conditions). 

When compared to the classical backgrounds the maxima in thermalized backgrounds 
are more pronounced, i.e. more localized (compare the scales in Fig. El ( c ) an d Fig. [HJ) • 
Accordingly, the IPR's are higher, varying between 1 = 4 and I = 11 for the lowest 
fundamental mode (cf. Fig. El (b)). Moreover, the regions outside the lumps have now 
a lower modulus. This can be quantified by the average of \<j)\, which for the thermalized 
background is lower than the classical one (and rises at the points of transition). For 
the lowest adjoint mode the IPR is even higher, I = 21 and I = 53 for periodic and 
antiperiodic, respectively, an effect which has already been observed in |17j . 

Thus, the overall behaviour of the lowest Laplacian eigenmodes in equilibrium back- 
ground fields is again analogous to the behaviour of fermionic zero modes [H] 5 . Still the 
lowest Laplacian mode seems to be broader than the fermionic counterpart, which is 
probably also the reason why the maxima are quite static, see Fig. 1161 

At this point it is desirable to discuss more quantitatively the feature of localization 
of the Laplacian modes. We have found it useful to add to the IPR measurements done 
so far a description suggested by Horvath 39 . It compares how quickly clusters (in our 
case of large modulus of the Laplacian mode) grow in size and how much - at the same 
time - they gain in norm, when their total volume increases by lowering the lower cutoff 
in the cluster definition. The size of the cluster has been defined here as the diagonal 
of the enclosing 3D cube (since the modes are almost static) and then divided by the 
diagonal of the 3D lattice. According to this description, a profile is called localized 
when small clusters already accumulate a large norm fraction. 

For the Laplacian modes we find mostly one big cluster. As Fig. El shows, the ad- 

5 Ref. .8, deals with SU(3) gauge fields (on a 20 3 -6 lattice) which should not differ in the general picture. 
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Figure 14: 'Hopping' of the modulus of the lowest Laplacian eigenmode. The rows rep- 
resent different modes. First row: periodic adjoint mode; second to fourth 
row: fundamental mode with £ = 0.020, 0.355 and £ = 0.480 representing the 
three (-intervals; fifth row: antiperiodic adjoint mode. The columns show 
different planes through the lattice. Left: plane through (22,24) = (12,1) 
with the maximum at (21,23) = (6,6) occurring in the first interval as seen 
in the second row; middle: plane through (22, 24) = (4, 4) with the maximum 
at (x±, 23) = (14, 14) occurring in the second interval as seen in the third row; 
right: plane through (22,24) = (6,3) with the maximum at (21,23) = (6,11) 
occurring in the third interval as seen in the fourth row; such that the plots 
contain the respective global maximum. The vertical scale is 0.1 for the 
adjoint and 0.05 for the fundamental plots, respectively. 
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Figure 15: Value and location (± one lattice spacing) of all local maxima that become 
global ones in some ^-intervals. 



joint modes are always local, but the fundamental modes are sometimes rather global, 
depending on the boundary condition £. This is the case for instance for £ = 0.300 (first 
transition point) and C = 0.355 (second (^-interval), cf. Fig. El (b) and (c). There the 
cluster has already its maximal size, i.e. has percolated through the lattice, although the 
accumulated norm is around three quarters only. 

Thus this analysis supports the view of Laplacian modes being close to wave-like (as 
has already been argued in the caloron background). Hence they should be able to reflect 
the background gauge field quite well everywhere on the lattice. This is an advantage 
for Laplacian gauge fixing and the filter method proposed in the next section. 

To end this section about Laplacian modes in equilibrium backgrounds, we want to 
mention a possibility to connect the behaviour of Laplacian modes in thermalized and 
classical backgrounds. It involves the smoothing of the gauge field background by ap- 
plying smearing (for instance 5 steps) to it. Then the maximum decreases in modulus 
(and can move or join with local maxima), whereas the minimum value increases. More 



Figure 16: Time dependence of the maxima in the three intervals, i.e. at £ = 0.020, 
C = 0.355 and ( = 0.480 with {x u x 2 ) = (6,12), (x 1 ,x 2 ) = (14,4) and 
(xi,x 2 ) = (6,6), from left to right, respectively. The vertical scale is 0.05 as 
in the corresponding space-space plots of Fig. El 
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Figure 17: Characterization of the cluster(s) occupied by Laplacian modes: accumulated 
norm (smooth curve) together with the relative size (step curve) as a function 
of the space-time filling fraction of the mode. The upper row shows (a) the 
lowest fundamental Laplacian mode with £ inside the first interval (the mode 
is localized), (b) at the first transition (where it has turned to global) and (c) 
inside the second interval (less global). The lower row shows the lowest adjoint 
Laplacian mode with (d) periodic boundary condition and (e) antiperiodic 
boundary condition (both are highly localized). 

importantly, the number of local minima decreases drastically and from some stage on it 
should be possible to use the minima as markers of the (gradually emerging) topological 
structure. However, such an approach mixes two techniques to suppress UV fluctua- 
tions and it would be cleaner to interpret directly the Laplacian modes of the original 
configuration. 

Interestingly, we have found an interrelation of the pinning of the Laplacian modes to 
a gluonic observable, namely the Polyakov loop. The maximum of the modulus of the 
periodic mode is correlated to a positive Polyakov loop, whereas the maximum of the 
antiperiodic mode prefers negative values, see Fig. 1181 (a). Of course, this expresses only 
an overall tendency, because there is a mismatch between the smoothness of the Lapla- 
cian mode and the roughness of the Polyakov loop. To make this statement quantitative 
one can measure the sum over the Polyakov loop weighted with \(f)\ 2 of the periodic and 
antiperiodic mode (the horizontal coordinate of the 'center of mass' of that figure). This 
gives weighted averages of 0.206 and —0.1154 for the two cases, respectively, while the 
ordinary Polyakov loop average is 0.005. 

We have smeared the gauge field which results also in smoothing the Polyakov loop, 
then we compared the latter to the Laplacian modes of the original configuration. The 
mentioned tendency becomes more pronounced, see Fig.^](b) and is now locally visible. 
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Figure 18: Scatterplots of the Polyakov loop (horizontal) at lattice sites where the modu- 
lus of the Laplacian modes (vertical) with £ = 0.02 (close to periodic, crosses) 
and £ = 0.48 (close to antiperiodic, circles) is large: (a) for the original con- 
figuration and (b) for the Polyakov loop measured on the configuration with 
5 smearing steps applied. 



As Fig. shows, the smeared Polyakov loop provides a collection of pinning centres 
for the modes to settle down at Polyakov loops near I2 and — 12, repectively. Hence 
the Laplacian modes on the unsmeared configuration know about the Polyakov loop 
landscape, which only emerges after smearing. 

The direction of this correlation agrees with the findings in classical solutions (compare 
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Figure 19: Comparison of contours of constant Polyakov loop (dashed curves) after 
smearing and of the unsmeared Laplacian modes at large modulus (70% and 
90% of the maximum, full curves), in lattice planes where the latter become 
maximal, (a) the periodic mode vs. positive Polyakov loop (0.7 and 0.9) in 
(^2,^4) = (12, 1), (b) the antiperiodic mode vs. negative Polyakov loop (—0.7 
and —0.9) in (x2, x&) = (6,3). 
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Fig. ^ (a) and Fig. |f)J) which seems to suggest calorons as 'underlying' the thermalized 
gauge fields. However, it can also be interpreted within the Anderson localization sce- 
nario, with the Polakov loop localizing the candidate 'minima of a random potential' to 
be occupied by the scalar or fermionic modes. 

5 A new filter method 

In the following we will introduce a new low-pass filter based on the Laplacian modes. 
It uses an exact representation of the links in the form of a sum of the latter, which, 
by truncation, should remove UV noise and find the 'underlying IR structures'. Hence, 
the spirit of our method is close to a Fourier decomposition or a (gauge invariant) high 
momentum cut-off. 

Technically, our approach is similar to the one in |18| . where the field strength was 
given in terms of a sum over fermionic modes. A similar mode truncation for the overlap- 
based topological charge has been used in [3U| and studied in more detail in |32] . 
However, we will obtain directly the link variables and in this 'reconstructed' configura- 
tion any observable can be measured. 

5.1 Derivation and properties 

We combine the definition of the gauge covariant Laplace operator, Eq. |T|). with a 
decomposition into its eigenmodes: 

M 

n=l 

At y = x + fi one immediately obtains 

<( x ) = ~ E WS.cfrXcC* + A) (12) 

n=l 

The corresponding formula for Up at y = x — fi is fully equivalent to Eq. (JT2|) , while 

1 N 

s ab = .E^A^ftW (13) 

n=l 

N 

= E A «.C<c( x Xc^) Vy^,x±/x (14) 

71=1 

could be used as a check for the approximation to be discussed now. Notice that the 
eigenfunctions enter these expressions in a combination that is invariant under a mul- 
tiplication of the eigenfunctions with a complex phase. More general, at a level with 
degeneracy k the expressions are invariant under a change of the basis (a global U(k) 
rotation acting on the index n). 
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The idea of the filter is to truncate the sum in Eq. hlfy) at a rather small number N 
of eigenmodes. When doing so, the question arises, how to relate the r.h.s. of Eq. (|12j) 
to a unitary link variable. We remind the reader that in the cooling method staples are 
added. This gives a quaternion, which becomes an element of SU(2) upon multiplying 
by a real number. We will use the same projection. However, the non-trivial task here 
will be to arrive at a quaternion in the first place. We will show now that the charge 
conjugation symmetry of the Laplacian helps to resolve this issue. 

For the <^>-terms in the sum we use the abbreviation 

<(x) ni? = < ? (x)< c (* + /i). (15) 

The charge conjugation symmetry © implies 

u'^(x) n _ c = e ^(x)* )f e -1 , e = ia 2 ■ (16) 

For Q = (periodic) and £ = 1/2 (antiperiodic) these contributions are automatically 
included in the sum in Eq. (|12j) (as every level is two-fold degenerate with eigenmodes 
(j) and (f)'). To make use of this relation for general C, we take the average of the link 
formula (|12|) over £ and — ( 



Z n=l 
1 N 

= -T. J2 X n,( [u^{x) ntC + eu^{x)* n ^e l ] . (17) 
1 n=l 

So far this manipulation is exact, but seems artificial. However, it will help in the 
truncation of the sum. Parametrising the matrix in the bracket we have 

: \ ) + ( s s ) - ( -:• ) ™ 

For a matrix of this form 6 the inverse is the same as the hermitean conjugate up to a 
factor, which is the determinant |u| 2 + |f | 2 and positive for all practical purposes. Hence, 
the matrix obtained by truncating the sum is an element of SU(2) up to the square root 
of the determinant by which we will divide. 

An alternative way to arrive at a unitary link is to use the unitary matrix in the 
singular value or polar decomposition. It agrees with our procedure under special cir- 
cumstances, for instance if the matrix to start with has a real determinant as is the case 
for C = an d C = 1/2- In our point of view the inclusion of both £ and — £ is very 
natural since these boundary conditions contain the same information of the Laplacian 
mode (same eigenvalue and charge conjugated eigenmode). 



J Eq. 1181 is actually the most general quaternion parametrised by Reu-l2+Imii-icr3+Re v-ia2+lmv 4ai. 
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To summarize, there is an exact formula for the link variables in terms of the Laplacian 
modes, namely Eq. Q17|). It is this sum we truncate to obtain the filtered links 



%(x)n,( = I -J2 Xn '( M^W + ^^c^ 1 ] ) M det=1 = M/VdetM , (19) 

V n=l / det=l 

where the operation (. . .)d e t=i forces the link variable to have unit determinant 7 (and 
we have dropped the factor 1/2). 

Eq. ()19|) is our final proposal for a low-pass filter acting on lattice configurations. It is 
a mapping from the original links U^{x) to the filtered ones U^{x)n,c, via the Laplacian 
eigenmodes with the phase £ in the boundary condition as a free parameter. The quality 
of the filter is controlled by N, where N = J\f reproduces the original configuration 
exactly (with no determinant correction necessary). Formally, the filtered links U^{x) 
look like composite fields, as they are produced by bilinears in (f>, Eq. (fT5|) . 

Our approximation keeps the full gauge covariance, because if a gauge transformation 
g{x) acts on (j)(x) from the left, then g^(x + /t) acts on cf>*(x + ft) from the right (and the 
determinant is gauge invariant) resulting in g(x)U /1 (x)g^ (x + fi) as it should be. 

It is instructive to investigate the extreme case of iV = 1, i.e. the contribution from 
the lowest mode only. Then the filter gives - independent of the starting configuration 
- a pure gauge U^(x), but with a Polyakov loop given by (. 

The best way to understand this is going into Laplacian gauge [El- That is to rotate 
the lowest mode to have only an upper nonvanishing component, which is real. Then the 
matrices Uu{x) would have only a real 11-component and adding the charge conjugate 
and normalising by the square root of the determinant makes all U^(x) equal to the 
identity. However, for the nontrivial boundary conditions of Eq. © one has to slightly 
rethink this procedure. The simplest way to define the gauge fixing then is to demand 
the form (real, 0) T for the (always periodic) ip^(x). By transforming back to (p((x), Eq. 
fj|. and plugging this into Eq. (|15[) the same considerations hold, with the only exception 
that u a ^ l (x) receives an additional factor exp(—2iriC / N^) . After charge conjugation 

and normalising the filtered links become U±{x) = exp(— 27ri^T^ / N^) . Hence the filtered 
Polyakov loop is 

-^V{x) N=u = cos(2vrC) Vx. (20) 

These considerations did not use that <j){x) is the ground state. Thus every eigenmode 
alone gives vanishing action and a Polyakov loop as above. Since the superposition in 
Eq. (j!9j) is nonlinear, we expect N = 2 to be the first nontrivial filter with nonvanishing 
action 8 . As we will see in the next section, the trace of the Polyakov loop will start to 
fluctuate in this case, but on average is still close to the value of Eq. f|2C7|) . 

7 We have used the positive root in Eq. 119H for all links. One can in principle also work with a local 
choice of the positive or the negative root, but we do not see a reason to introduce a Z2-freedom 
here. 

8 This is true unless C, — or £ = 1/2 for which the two- fold degeneracy makes N = 3 (equivalent to 
iV = 4) the first nontrivial case. 
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Let us add a few remarks on the possible stability of the filtered links U u (x)n^ under 
variations of N. From the normalization of (f> n (x) follows that the entries of u n (x) are 
roughly 0(1/ Af). The prefactor X n is rising slightly with n, such that one can expect 
the terms in the sum to be of the same order of magnitude. 

The only exception to this argument is the caloron background at £ = 1/4, where 
the lowest Laplacian eigenvalue is strongly suppressed w.r.t. the other eigenvalues. As 
a result, the lowest mode practically does not contribute to the filtered links. Here 
the situation is analogous to the fermionic filter in ^H], where the zero mode does not 
contribute. 

Actually, the staggered symmetry of the Laplacian can be used to improve the con- 
vergence. Relation © between low and high modes immediately gives 

u "n( x )N-n,c = -Un(x) n ,<: ■ (21) 

Therefore, the inclusion of the upper end of the spectrum does not change the contribu- 
tions to the filtered links locally. It only reweights the terms in the sum ()19[) 

- X n -> AD - 2X n . (22) 

With this insertion one has to sum in Eq. (|19|) only over half the spectrum to obtain the 
original link. Now each subsequent term u u (x) n has a smaller weight. In particular, the 
weight of the ground state is the biggest. 

The price to pay is to include the upper end of the spectrum of the Laplacian, which 
a priori contradicts the meaning of a low-pass filter. We decided not to include the 
highest modes and to stick to Eq. (|19|h Nevertheless, we have checked the consequences 
of such a modification for the results presented below. The changes are quite small 
apart from the case of reconstructing the caloron background from N = 2 modes at 
C = 1/4, the case discussed above. Therefore, it seems that the local information u^(x) n 
is more important for the filtered links than the relative weights A n of the eigenfunctions. 
One might speculate whether eigenmodes other than the lowest ones could be used and 
whether the weights could be chosen arbitrarily (constant or random). 

5.2 Classical objects seen through the filter 

In this subsection we test the filter in a controlled environment, namely for calorons as 
smooth configurations. 

First we want to study a technical detail, namely how big the determinant is that 
appears in (j!9[) to scale up the sum to an SU (2)-valued link. In other words, how much 
of the link is captured by the superposition of a finite number of modes (in a naive sense 
without inspecting any structure in the links). 

From Tab. one can read off that this determinant is small (the correction factor is 
large) when only a few modes are taken. This is to be expected from the smallness of 
the 4> n (x) that are square-normalized on the entire lattice. The determinant grows with 
N. Of more relevance is that the minimum and the maximum of the determinant taken 
over the lattice approach each other, such that all links appear 'equally well filtered'. 
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N 



1 4 10 50 200 



miiLj. log 10 det -15.3 -12.1 -8.5 -6.1 -4.1 
maxj; log 10 det —11.9 —8.6 —7.6 —5.5 —3.9 

Table 1: Behaviour of the determinant used to project to an SC/(2)-element in the trun- 
cated sum (|19|). depending on the number N of modes. The background is the 
large caloron on 16 3 • 4 with £ = 1/4. 

The most important question in the caloron context is whether the filtered links pro- 
duce the action density lumps of the two constituent monopoles including the typical 
structure in the Polyakov loop. In Fig.|2niwe give the observables corresponding to Fig. 
n (a) in the filtered configuration with different £ and N. 

For the chosen N -C M, the boundary condition parameter governs the average 
Polyakov loop. In the intermediate case £ = 1/4 the average Polyakov loop trace is 
with extrema of nearly ±1 at the monopole cores, see Fig. 1201 (b). In the periodic 
case C = the Polyakov loop is almost everywhere I2 with a dip resembling V = — 12 
at the corresponding monopole, see Fig. 1201(a). The antiperiodic case (not shown) is 
complementary, V ~ — 12 with a signal at the other monopole. 




(a) (b) (c) 

Figure 20: Action density, topological density and Polyakov loop computed from filtered 
links in the background of the large caloron of Fig. ^ (b) ( = 1/4 with 
N = 4 modes (action and topological density multiplied by 100); (c) more 
modes N = 150 (scale 400 just as in the original plot Fig. H](a)). (a) periodic 
boundary condition £ = with N = 4 modes (scale 30). 

This picture stays the same when more modes are taken into account 9 . The Polyakov 
loop for the periodic case develops a stronger dip then, but is still very close to I2 
anywhere else. For the intermediate case the average stays close to 0. Locally it agrees 
almost perfectly with the one of the original configuration, see Fig. |20] (c) for N = 150 
vs. Fig. ^(a). To make this more quantitative we give the values of the Polyakov loop 
at the action density lumps (see below) and on average in Tab. |2J 

The findings for the reconstructed action and topological charge are analogous. In or- 
der to compare the two, we have computed them with the help of an 0(a 4 ) improved field 



Although, of course, in the limit N — > M the Polyakov loop is the original one for all (,. 
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1 n 

1U 


1 nn 

1UU 


1 


zuu 


original 


X3 at lumps 


4.0 


3.7 


4.1 


4.3 


3.7 


4.4 




13.0 


13.3 


12.9 


12.8 


13.3 


12.6 


action density 


0.014 


0.010 


0.0023 


0.0037 


0.002b 


0.0032 


top. density 


0.013 


0.010 


0.0019 


0.0035 


0.0026 


0.0032 


Pol. loop 


±0.84 


±0.80 


±1.00 


±0.99 


±0.99 


±0.97 


total action 


2.4 


2.2 


1.62 


1.34 


1.63 


1.07 


total top. charge 


0.98 


0.99 


1.00 


1.00 


1.00 


1.00 


average Pol. loop 


-2- 10~ 6 


-4- 10~ 6 


-0.02 


0.0001 


-0.0002 


-0.0003 



Table 2: Quality of the reconstruction of a large caloron from the filter with different 
number of modes. Both, local quantities (interpolated) and global quantities 
converge to the original values within some error margin. 

strength tensor |43| . The periodic case sees only one monopole 10 , while the intermediate 
case is able to detect both (Fig. |2UJ). We conclude that the boundary condition ( = 1/4 
is the appropriate one when one wants to describe calorons with maximally nontrivial 
holonomy through the filter, as is also clear from the discussion of the contribution of a 
single mode in the last section. 

Remarkably, the filter with the number of modes as low as ./V = 4 already has quite 
some knowledge about the classical structures. Although the action density lumps are 
a bit spiky and their locations are not perfect as recorded in Tab. [2] (notice also the 
different scales for the action and topological density used in Fig. [30] (b) vs. (c) and 
Fig.n(a))) they clearly reflect the constituent structure with opposite Polyakov loops. 
The two lumps are almost selfdual (cf. third and fourth row in Tab.|2J). Furthermore, 
the adopted definition of the topological charge works quite well for the filtered links 
and gives a total topological charge close to 1. We stress that the filtered configurations 
have not undergone further cooling and their approximate selfduality is a remnant of 
the original configuration passed on by its lowest-lying Laplacian modes. 

The inaccuracies of the N = 4 case are cured by increasing the number N of modes 
used in the filter, however, not systematically. The quite perfect case N = 150 used 
for the plot in Fig. [20] ( c ) is contrasted by N = 200, where some of the observables 
are further away from the original configuration. We interpret this as a limitation on 
the precision of the filter. Away from the ultimately exact limit N = M, the filter will 
reproduce the classical background only within some error and hence it is not useful 
to take more than a few hundred modes into account (at least for this example). This 
does not spoil the use of the filter at all, since its intended application will be mainly to 
thermalized configurations (see below). There it is not the aim to reproduce the rough 
background to a high precision, but to keep only a minimal structure of it such that it 
still captures the physical randomness. 

We note in passing that the Taubes winding inside caloron configurations as well as 
the ring structure for the charge 2 example are displayed by the filtered configurations, 
too. 

10 Increasing N makes the other monopole visible, but with a height much lower than that of the first 
monopole. Thus this choice of £ reproduces the equal mass constituents quite badly. 
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5.3 Are the filtered fields confining? 

The next two subsections are devoted to the application of the filter to thermalized 
configurations. One of the main issues here is whether the filtered links still give rise to 
confinement. 

To this end we have measured the Polyakov loop correlator on the ensemble of 50 con- 
figurations on a 16 3 -4 lattice created at /3 = 2.2. The Polyakov loop as the deconfinement 
order parameter has average and therefore we choose £ = 1/4 in the construction of 
the filter, Eq. (|19|) . Indeed, the filtered Polyakov loop (tr'P/2) averaged over the lat- 
tice has an expectation value compatible with (with 50 configurations giving standard 
devitations from 0.015 for N = 2 to 0.083 for N = 100 whereas the original standard 
deviation is 0.025). This again confirms that the average Polyakov loop for small N 
follows the one of the trivial case N = 1. In a way, we adapt the Polyakov loop average 
by a parameter in our filter and then look at correlations in fluctations on top of it. 

Fig.r^Jshows how the Polyakov loop looks locally (for a fixed configuration and lattice 
plane). When taking more modes into account, the Polyakov loop deviates further and 
more often from the average 0. Fig. FJ3 makes this statement more quantitative. It 




(c) (d) 

Figure 21: Polyakov loop 'evolution' with increasing number N of modes for a thermal- 
ized configuration in a fixed lattice plane (the one of Fig. I25|) . (a) N = 2, (b) 
N = 10, (c) ./V = 100 and (d) the original configuration. The vertical axes 
are from —1 to 1. 
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displays the distribution of Polyakov loops for one filtered configuration . As is clear 
from the discussion so far, the distribution for N = 2 is quite narrow and broadens with 
growing N. At N = 100 it is very close to the one of the original distribution and the 
Haar measure \J\ — (trV/2) 2 . 

0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


-1 1 

Figure 22: Distribution of Polyakov loops over the lattice sites for one configuration 
when filtered with the intermediate boundary condition £ = 1/4 for N = 2 
(narrow distribution around 0, divided by 15 to fit in the same plot) and 
N = 100 (solid line) compared to the orignal configuration (dashed) and the 
Haar measure (dotted). 

Fig. 1231 shows the logarithm of the Polyakov loop correlator, related to the interquark 
potential, measured at the filtered ensembles with N = 2 to N = 100 together with the 
unfiltered one. There is clear evidence that the filtered Polyakov loop correlator decays 
exponentially with the distance, in roughly the same window as the original one does. 
The potential after filtering has no sign of a Coulomb regime since the filter has washed 
out short range fluctuations as it should. The filtered curves are also shifted vertically. 
The value at zero distance represents the width of the Polyakov loop distribution. It 
approaches the original one from below, because the filtered distributions are narrower, 
as was shown in Fig. 1221 

In order to cast the confining behaviour into numbers we have performed an exponen- 
tial fit c exp(— or) in the range between 3 and 6 lattice spacings 12 . The slope a is directly 
proportional to the string tension a = a/N^a 2 . We will use the corresponding string 
tension of the original configurations 13 as the reference observable. The lower end of the 
fit region was taken such that the numerical value of a for the original configuration is 
stable when compared with a fit over 1 to 6 lattice spacings including a Coulomb part. 
Towards large distances the fit was limited by statistical errors. The obtained values for 

lx We thank A. Wipf for suggesting to show the ^-dependence of the local Polyakov loop distribution. 
12 Our lattice spacing is 0.21 fm provided cr(0) = (AAOMeV) 2 . 

13 The obtained value of the string tension (for the original configurations, see Tab.|^J is equivalent to 
0.78 er(0) whereas a parametrization a = ^Jl — (T/T c ) 2 <r(0) 1111 would predict 0.66 <r(0). 
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Figure 23: Polyakov loop correlator plotted on an inverse logarithmic scale over the 
distance for different N compared to the original one. 

the string tension and estimates of its error are given in Tab. Although the approach 
of the filtered string tensions for the number of modes between N = 2 and 100 to the 
original one is again not monotonic, the latter is reproduced within 15% (for N = 4 
modes the reproduced string tension is almost perfect). 

A priori it is not obvious whether the Polyakov loop fluctuations in the filtered and 
the original configurations are correlated or whether the filter creates independent fluc- 
tuations (a pointwise scatterplot of the two respective Polyakov loops against each other 
shows a correlation setting in not below N = 50 modes). A similar question has been 
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Figure 24: Test whether the original and filtered configurations contain correlated fluc- 
tuations disordering the Polyakov loop. Correlators C(r) (upper curve) and 
Cind( r ) (lower curve), see text. 
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investigated in [15] concerning two Wilson loops W\(C) and W2(C) 14 . There the corre- 
lation {Wx{C)W%{C)) (i.e. between different Wilson loops along the same curve C) was 
compared to {Wi(C)) (W2(C)) representing the hypothetic independence of the Wilson 
loops (and decaying exponentially with the area of C). Adopting this idea to our case 
we have to compare 

C(\x-y\) = ^-(trV(x)trT(y)-tTT N (x)trT N (y)) to (23) 
lb 

C ind (\S-y\) = ^-(trV(x)tvVm(^V N (x)trV N (y)). (24) 
lb 

We did this for N = 4 and show in Fig. |21|that the correlator C(r) becomes independent 
of the distance r rather quickly. Ci n d(r) of course decays exponentially with the sum of 
the string tensions of the original and the filtered case. We conclude that the fluctuations 
disordering the Polyakov loop for the original and the filtered configuration are not 
independent but rather correlated. 

5.4 Structures found by the filter 

Finally, we give a first description of the vacuum structures that appear when the filter is 
applied to an individual configuration from an equilibrium ensemble, representing finite 
temperature in our case. 



N 


2 


4 


10 


50 


100 


original 


a/10 5 MeV' 2 


1.732(31) 


1.447(25) 


1.299(17) 


1.385(13) 


1.574(25) 


1.509(70) 


Wilson action 


6b 


84 


10b 


177 


234 


2131 


max. # clusters 


37 


50 


70 


79 


104 


342 



Table 3: Observables to characterise the filter when acting on equilibrium configura- 
tions. The string tension is computed from 50 configurations, the other two 
observables are for the example configuration characterised locally in Fig. 1211 
and ESI 



The behaviour of the Polyakov loop was shown in the previous section (Fig. I21|). The 
total action of the filtered configuration behaves like iV 1 / 3 in a range from N = 2 to 
N = 100 eigenmodes. The interpretation of this particular dependence is not obvious, 
but the overall picture is clear: the more modes are included in the filter, the more 
fluctuations occur and contribute to the total action. In the example we use (cf. Fig. 
1211 and l25|) the action in the first nontrivial case N = 2 is b6 in instanton units of 2ir 2 , 
compared to a total plaquette action of 2131 for the original configuration. The first 
number should be read as the minimal content that survives the filter. The only way 
to achieve an 'even smoother configuration' would be to apply the filter to the already 
filtered configuration. 

The global maximum of the action density is only weakly varying with N, meaning 
that there are peaks of about the same height at every stage of the filter. From Fig. 1251 

14 We thank J. Greensite for urging us to perform the following test. 
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Figure 25: Left column: 'Evolution' of the action density with growing number of modes 
in the filter, (a) N = 4, (c) N = 10, (e) N = 50 and (g) N = 100. Chosen 
is the lattice plane (x2,x^) = (6,4) which contains the global maximum for 
N = 10 sitting at (x\,xz) = (10, 13). For comparison the action densities for 
the same configuration after 5 steps of smearing (b), after 2 steps of cooling 
(f) and for the original configuration (h) are presented with uniform scale. In 
the second row a comparison is made between the action density (c) and the 
topological density (d), both after filtering with N = 10. 
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showing the action density in a fixed lattice plane for different N, it is evident that these 
peaks are quite narrow. This effect is counterintuitive, but has also been observed for 
caloron backgrounds in the last section. The peaks look the same in a space-time plane, 
too, i.e. they are not static. The figure also shows that at least some of the peaks are 
stable w.r.t. N, that is they stay as local maxima. Their width is not changing much 
either. 

In order to more quantitatively describe the filtered structure we have performed a 
cluster analysis. Lowering the threshold we have recorded the number of clusters and 
their respective volume, size and accumulated action. The number of clusters first raises 
and after reaching a plateau decreases again due to cluster mergings, while the peaks 
remain. 

As shown in Fig. 1261 the corresponding curves fall almost on top of each other, when 
the number of clusters is taken relative to its maximum at that N (see Tab. El . That 
suggests that the difference in the filtered action for different N is due to the different 
number of clusters. Otherwise they are rather similar. Other indications for that are 
the volume, the size and the action per cluster, which - plotted as a function of the 
threshold - fall on top of each other as well. 

In Tab. El we give for each iV the maximal number of clusters. The latter appears at 
a threshold of roughly a third of the global maximum, where all clusters together have 
accumulated a few percent of the lattice volume and roughly 15% of the action. It is 
interesting to notice that the maximal number of action density clusters roughly agrees 
with half the total action (in instanton units, see Tab.Elthird and second row). A lower 
bound for the estimated action per peak is thus 0.3. More realistic is to assign part of 
the remaining action below the threshold to the peaks (sitting in their tails). From Fig. 

we read off that at a threshold of 0.01 (around 1/15 of the global maximum) the 
number of clusters starts to fluctuate because the region of action density noise has been 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 



Figure 26: Relative number of action density clusters as a function of the lower threshold 
for a filtered configuration at different N. 



35 



reached. At that threshold the accumulated action of the clusters is roughly 60% of the 
full one. This results in an average action per peak of 1.2 in instanton units. Hence, 
from an action density point of view, this estimate is compatible with an interpretation 
of the peaks as (possibly fractional) instanton lumps. 

Indeed, the local maxima of the action density are often accompanied by local extrema 
in the topological charge density. This is visualized for N = 10 in Fig. 1251 fd). However, 
the topological charge density does not equal the action density at the peaks. For the 
global maximum at N = 10 the ratio of both quantities is 0.87, but for other local 
maxima no signature of (anti)selfduality was found. Moreover, the total topological 
charge of generic configurations expressed in terms of the filtered links is not close to 
an integer. In this respect, the configuration, although it is filtered, is not smooth 
enough to make the 0(a ) improved field strength definition of the topological charge 
work. Fermionic definitions via the index should be used to determine whether the filter 
preserves the total topological charge. 

The peaks seem to have nothing to do with the maxima of the modulus of the Laplacian 
mode in that background, but with the determinant used as the normalization factor in 
Eq. (|19|) . Obviously this determinant is given in terms of eigenmodes (and eigenvalues) 
of the Laplacian, too, thus it also contains information about the gauge background. 

Another interesting question is whether the new filter is related to smearing and 
cooling. In Fig. we compare the corresponding action densities to the filtered ones, 
where the number N of modes was chosen such that the total actions are comparable. For 
instance, 2 cooling steps result in an action of 172, which therefore compares to iV = 50 
(see Fig. |3] second row) and 5 smearing steps result in 84 which compares to N = 4. 
As the third and the first row of Fig. 123 show, the corresponding peaks seem to agree 
locally, which is a nontrivial correlation, because both methods lower the total action 
by a smoothing procedure, but in completely different ways. In contrast, 5 cooling steps 
lead to an action of 56, which seems to be comparable to N = 2. The action density, 
however, looks very different. It is clear that for more cooling steps the correlation to 
the filter has to break down, since cooling typically drives towards action plateaus of a 
few instanton units and finally removes the string tension. This is avoided in the filter 
method. 

More work has to be done to better understand the (iV-dependent) spiky structures 
induced by the filter. 

6 Discussion and outlook 

We have investigated Laplacian eigenmodes in both classical and thermalized gauge 
configurations with respect to their capability to analyze certain properties of the back- 
ground. Our localization studies for classical backgrounds lead to the conclusion that 
there is an analogy to fermionic (near zero) modes. In caloron backgrounds they de- 
tect the monopoles by a minimum and a maximum and hop between these constituents 
upon changing the boundary condition angle £. We conclude that such a localization 
is merely due to information about the gauge background in the covariant derivative, 
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which is modified by the angle £ in the same way, and does not depend much on the 
spin of the analyzing field. 

Yet, there is some difference between fermionic and Laplacian modes. This concerns 
for instance the behaviour under intermediate boundary conditions. For the caloron 
fermionic zero mode one lump grows at the expense of the other, while the Laplacian 
ground state develops a valley. Also the precision of the localization is lower than for 
the fermions: the maximum of the modulus is less pronounced and static even for a 
time-dependent action density of the background. For the minimum these problems 
do not occur, but the locations of both, minimum and maximum, deviate from the 
constituents. Moreover, the modulus of the mode approaches the average value away 
from those structures and hence the IPR's are rather small. That is why the Laplacian 
modes resemble modified waves rather than exponentially localized discrete states (in 
continuum language). 

The wave character also applies to the Laplacian modes in the adjoint representation, 
where the underlying structures (constituent monopoles or instantons) are detected by 
minima only 15 . We expect them to become zeros in the continuum, which is natural from 
the point of view of the Laplacian Abelian projection. For antiperiodic adjoint modes 
these minima form even two-dimensional sheets between the constituent monopoles. 
We have illustrated this for a semiclassical background, too. Signatures of the Taubes 
winding can be found in both representations. 

In thermalized backgrounds, there is clear evidence that the lowest mode changes 
its global maximum with the boundary conditions. Different intervals in the boundary 
condition angle £ emerge, where different local maxima take over the role of the global 
maximum. We have observed up to 4 jumps per configuration and the corresponding 
lattice locations seem to be close to randomly distributed. These locations are also 
visited by excited Laplacian modes and those in the adjoint representation. 

The corresponding global minima are not stable under the boundary condition (and 
the number of local minima is large). Therefore they could not be used as a practical 
tool for localization. Thus, a straightforward interpretation of Laplacian modes on 
thermalized backgrounds in terms of classical objects is not possible. In this context it 
would be interesting to study the effect of quantum fluctuations on the Laplacian modes 
by heating a classical background. 

Actually, the IPR is basically insensitive to minima. Another observation that ques- 
tions the use of the IPR is that, in most cases we studied, it is proportional to the value 
at the global maximum. This large value apparently dominates the sum in the IPR 
definition and therefore the latter gives no new information. In addition to the IPR 
measurements we have performed a cluster investigation. Somewhat unexpectedly we 
have found that the fundamental Laplacian mode can even be characterised as a global 
structure (depending on £). 

Apart from the weaker localization and their more static nature, the Laplacian modes 
behave again similar to fermionic modes. A natural next step is to clarify whether 
Laplacian and fermionic modes on the same configuration see the same locations. In 

15 For adjoint fermions there are more zero modes, with maxima at the constituents 0- 



37 



order to eventually go beyond a purely empirical description it is desirable to sort out 
and measure the relevant gluonic features in the gauge field backgrounds, for instance 
the topological charge density. 

Smoothing techniques like smearing could be used for this purpose, although they have 
the disadvantage to modify the gauge background. Then the Laplacian modes become 
similar to those in classical backgrounds: they level off thereby lowering (and joining or 
slightly moving) the lumps and smoothing the minima. 

We have pointed out an interesting relation of the unsmeared Laplacian mode to a 
gluonic observable after smearing: the smeared Polyakov loop provides pinning centres, 
some of which the Laplacian mode occupies. Such a finding in itself is in accordance 
with both hypotheses, calorons underlying the thermalized configuration and Anderson 
localization. More work has to be done to better assess the real mechanism. 

Another interesting question is whether the Laplacian modes reveal signatures of the 
deconfinement phase transition. As a first signal we have seen a qualitative difference 
in the C-dependence of the spectrum. Of course the Laplacian modes should also be 
investigated at zero temperature, to see to what extent the properties described in this 
paper remain. 

Laplacian eigenmodes are the natural ingredients for a Fourier-like filter. We have 
described a new method to reconstruct the link variables by truncating a sum over 
Laplacian modes. Because of the use of eigenmodes one might view this technique as a 
nonlocal smoothing. It is actually not too expensive, the computation of N = 10 modes 
on a 16 3 • 4 lattice including the reconstruction of the links takes a minute on a 1.7 
GHz PC. It should be stressed that the filtered links allow for the measurement of any 
quantity and that the procedure is not biased towards any particular degree of freedom 
in the QCD vacuum. 

The striking properties of the filter are that it reproduces classical structures, in 
particular selfduality, and preserves the string tension (within 15%) when applied to 
equilibrium configurations. The number of Laplacian modes can be kept remarkably 
low, typically 16 4 out of more than 10 5 eigenmodes start to reproduce the mentioned 
features qualitatively. Because the filter keeps the relevant long-range disorder, we claim 
a 'low mode dominance' in the confining properties of lattice gauge theory. 

The filter has a lower limit, namely N = 2 modes. It might well be that the content 
of action density and Polyakov loop fluctuations for this case is minimally required to 
keep the long-range physics. For the same reason it is clear that the observed similarity 
of the filter with smearing and cooling is limited, for instance to early stages of cooling. 

Apart from the number of modes the only filter parameter is the phase in the boundary 
condition. We have shown that it influences the quality of the filter and that in the 
confined phase (at finite temperature) it is best set to the intermediate value £ = 1/4. 

Concerning the tomography of the filtered configurations we have not yet reached 
a final understanding. The filtered Polyakov loop is narrowly distributed around 



This number, referring to (3 = 2.2 on a 16 3 • 4 lattice, might depend on the lattice volume and the 
coupling constant. 
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and approaches the original distribution for iV ~ 100. In the emerging action density 
isolated peaks appear, both for caloron and equilibrium configurations. This seems 
counterintuitive, taking into account that the filter uses the lowest Laplacian modes. 
One should keep in mind, though, that the contribution of the high end of the spectrum 
is almost the same as for the low end (due to the staggered symmetry). 

It seems that the peaks are correlated to the normalization factor (inverse square root 
of the determinant) that has to be applied to project the filtered link back to SU(2). 
This mechanism as well as a possible relation to the structures found by Fourier-filtering 
in Landau gauge [E] and to singular gauge fields jUj has to be clarified. 

In this context it might be helpful to improve the filter procedure. First of all, the lat- 
tice Laplace operator could be replaced by an improved version. The second opportunity 
is to interpolate the eigenmodes. This leads to filtered links on finer lattices (similar to 
inverse blocking |48[ I49|). which can be subject to blocking. A prerequisite for this idea 
is that not only the modulus but also the components of the lowest-lying eigenmodes 
are fairly smooth. This is actually the case in Laplacian gauge. 

There are some possibilities to generalize our method or to apply it in another physical 
context. The obvious applications are to zero temperature and to the deconfined phase, 
respectively. For the latter we have checked that filtering with £ = 1/4 results in no 
string tension, however, one might be forced to fix £ in a different way. The quality of 
the filter in this phase could be checked by virtue of the spatial string tension. 

The generalization to higher gauge groups is not completely trivial. We have used 
the charge conjugation symmetry (relying on the pseudoreal nature of SU(2)) in the 
derivation of the filtered links becoming unitary. For gauge groups SU(N) the singular 
value decomposition seems to be the only alternative. 

Fermionic modes could also be used to reconstruct the gauge field, provided one is able 
to project out the spin indices to arrive at an exact formula for the links. This variant 
of the filter will be more expensive, but might be advantageous concerning topological 
properties of the gauge field background. 
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